Nonlocal effect on resonant radiation force exerted on semiconductor coupled quantum well nanostructures
Zhang Jin-Ke1, 2, Zhang Ting-Ting1, 2, Zhang Yu-Liang1, 2, Wang Guang-Hui1, 2, †, Deng Dong-Mei1, 2
Guangzhou Key Laboratory for Special Fiber Photonic Devices, South China Normal University, Guangzhou 510006, China
Guangdong Provincial Key Laboratory of Nanophotonic Functional Materials andDevices, South China Normal University, Guangzhou 510006, China

 

† Corresponding author. E-mail: wanggh@scnu.edu.en

Abstract

Based on the microscopic nonlocal optical response theory, the resonant radiation force exerted on a semiconductor-coupled quantum well nanostructure (CQWN), induced by the nonlocal interaction between lasers and electrons in conduction bands, is investigated for two different polarized states. The numerical results show that the spatial nonlocality of optical response can cause a radiation shift (blue-shift) for the spectrum of the resonant radiation force, which is dependent on the CQWN width ratio, the barrier height, and polarized states sensitively. It is also confirmed that the resonant radiation force is steerable by the incident and polarized directions of incident light. This work may provide an advantageous method for detecting internal quantum properties of nanostructures, and open novel and raising possibilities for optical manipulation of nano-objects using laser-induced radiation force.

1. Introduction

Utilizing optical force to alter the motion of micrometre-sized particles[1] and neutral atoms[2] could have potential applications in the manipulation of microscopic particles and of individual atom. Therefore, the optical force and optical manipulation of nano-objects have been attracting extensive interests and have been successfully applied in various fields such as plasmonic nanoparticles, quantum dots, biology, etc.[311] For example, Nieto-Vesperinas et al.[5] had announced the study of the optical force on a small particle with both electric and magnetic response, and clarified the origin and significance of the electric–magnetic dipolar interaction force, which is related with the angular distribution of the extinction cross section and of scattered light.

Recently, the nonlocal optical effects in semiconductor and metal nanostructures have been drawing more and more attention.[1220] The so-called nonlocal optical response is that the polarization at a spatial point is induced by the applied optical fields not only at the same point, but also at other positions within the extent of the relevant wave function of electrons. Therefore, when we take the nonlocal optical response into account, susceptibility tensors and dielectric permittivity would be more complex. These pioneering works of spatial nonlocality in semiconductor materials have developed into two successful research aspects: nanostructures[21] and metamaterials.[22] For instance, Mortensen[23] had found some important nonlocal effects in metamaterials and Batke et al.[24] demonstrated a strong nonlocal interaction of two-dimensional (2D) magnetoplasmon resonance with harmonics of the cyclotron resonance in AlGaAs/GaAs heterostructures. Moreover, there are also a great deal of interests in nonlocal optical effects for mesoscopic systems, resulting in the development of nonlocal theory and nonlocal optical properties in plasmon and semiconductor nanostructures.[2528] The nonlocal optical response of metal nanostructures with arbitrary shape had been investigated and the particular importance of nonlocal effects in apex structures had been indicated by McMahon et al.[20] Furthermore, Kosionis et al.[29] had also taken nonlocal effects into consideration in the research on optical effects in quantum-dot/metal–nanoparticle hybrid systems. As far as we know, the spatial nonlocal effects on optical force have not been clarified. It is reasonable to believe that the spatial nonlocal effects may have an important influence on the optical force exerted on nanodevices. The physical motive of our paper is to demonstrate and explain the optical force exerted on a coupled semiconductor quantum well nanostructure (CQWN) from the resonant intersubband transitions for two different polarized states under spatial nonlocal effects.

This paper is organized as follows. In Section 2, a basic theoretical framework for optical force applied on a CQWN with two kinds of polarized states is presented by taking the spatial nonlocality of optical response into account. Section 3 is dedicated to the results and discussions for AlGaAs/GaAs CQWN. And finally, our conclusions are given in Section 4.

2. Theoretical framework for optical force on a CQWN with nonlocal response

In what follows, let us consider that a monochromatic plane wave with an incident angle θ and an angular frequency ω in the cases of p- and s-polarized states, is incident on a CQWN confined in the z direction, as shown in Fig. 1. In general, the Lorentz force acting on an object induced by an electromagnetic field is described as[30]

where and are the electric and magnetic fields, respectively. is the charge density, and is the current density. V is the volume of the force body. The charge density can be written as , where is the induced polarization intensity.

Fig. 1. (a) The schematic diagrams of radiation forces acting on a CQWN confined in the z direction, where Fx and Fz denote the x and z components of the radiation force induced by a monochromatic plane wave with p- or s-polarized state, respectively. (b) The schematic diagrams of the CQWN, where LL, LB, and LR denote the left well, mid-barrier, and right well widths, V1 and V2 denote the mid-barrier and outer barrier height, respectively. The signs and denote polarization direction of p- and s-polarized light waves, respectively.

For time-harmonic electromagnetic fields, we have

where , , and are complex functions of position and frequency. Re denotes the real part. The time-averaged force in the case of time-harmonic electromagnetic fields is derived as[8]
where * denotes the complex conjugate. Under the nonlocal linear response, the resonant part of the induced polarization is described as
where is the nonlocal susceptibility tensor. Due to the system with translational invariance in the x direction and confinement in the z direction, the complex electric field and polarization can be assumed to have the forms and , where kx stands for the wave number in the x direction.

In the following, let us consider the AlGaAs/GaAs CQWN restricted in the z direction, as shown in Fig. 1. Based on the effective mass approximation, the Schrödinger equation for electrons in CQWN can be written as

with the confinement potential V(z)
where z represents the confinement direction of the CQWN; V1 and V2 denote the mid-barrier and outer barrier height, respectively; is Planck constant; is the effective mass of an electron in the conduction band; and denote the eigenenergies and eigenfunctions of the CQWN system, respectively. Due to the translational invariance in the xy plane and the confinement in the z direction of the system, and can be written as
and
here, and are electronic wave vector and position vector in the xy plane. The solutions of one-dimensional Schrödinger equation , with , are and , which are named the envelope wave functions and the corresponding transverse energies of the n-th subband, respectively.

For the finite-depth CQWN, the transverse energy eigenvalue and the corresponding electronic envelope wave function in the CQWN can be expressed as

where , , and . We can obtained the coefficients an, , , , , , , en, and the transverse energy eigenvalue by using the continuity of envelope wave function at the boundaries.

By substituting Eq. (3) into Maxwell’s equations, we can derive the integral–differential equation for the self-consistent complex electric field under the nonlocal linear optical response,[15] viz.

where the tensorial operator
Here c is the speed of light in vacuum. denotes the unit tensor. and denote unit vectors in the x and z directions, respectively. is the relative dielectric constant of the background medium. By the one-body density matrix method,[15] the components of the linear polarization, , for the two-level CQWN model can be obtained in the low-temperature ( ) and long wavelength approximation ( ) limits as follows:
where
and , e is the electron charge; is the off-diagonal relaxation rate; is the Fermi energy of the CQWN system, where Ns is the donor surface concentration.

By the Green’s function method,[15,21] one can solve the components of the microscopic electric field ( for s-polarized state, and for p-polarized state) from Eq. (9) as follows:

where is the incident field, and is the wave number in the z direction. In addition, , and , where . are the elements of the tensorial Green’s function of the CQWN system[15]
where sgn(z) is the usual sign function.

By multiplying both sides of Eqs. (12)–(14) by and , respectively, then integrating the three equations over z across the CQWN, one can obtain

where

Note that, for the following analyses, equations (16)–(18) can be written in vectorial form by the rotating wave approximation[21]

where is a 3 × 3 diagonal matrix with the elements given below[21]
The vectors and , where , , and .

However, if we only take local optical response into account, the complex electric field and polarization can be written as

where kx and kz stand for the wave number in the x and z direction in the CQWN, respectively; is the local linear susceptibility tensor, whose elements , , here and L is the total width of the CQWN.

By substituting Eq. (20) and Eq. (21) into Eq. (2), we can obtain the numerical solution of the resonant radiation force exerted on the CQWN in the case of the local optical response.

3. Results and discussion

To discuss the optical force exerted on the CQWN, we take the spatial nonlocality of optical response into account in detail. Some parameters used in the following calculation are adopted as:[12,14] (m0 is the mass of a free electron), Ns=1011 cm−2, , and the mid-barrier width LB=1 nm. The incident intensity is assumed to be 0.05 W/cm2.

In Section 2, based on Eqs. (10)–(14), we can solve the components of the microscopic electric field and induced polarization by the Green’s function method[15,21] and one-body density matrix method,[15] respectively. Substituting these components into Eq. (2), we can obtain the numerical solution of the resonant radiation force, then analyze the resonant radiation force spectra.

In Figs. 2(a) and 2(b), the resonant radiation force is shown as a function of the normalized photon energy for asymmetric CQWN of two different width ratios between the right well width LR and the left well width LL: LR/LL=2.5/5 (black line), LR/LL=7.5/5 (red line) with the left well width LL=5 nm in two cases: (a) p-polarization and (b) s-polarization with the incident angle , and the barrier height V1=V2=2.0 meV. For the sake of contrast, we also plot the resonant radiation force of the symmetric CQWN for three different widths: LR=LL=3 nm (black lines), LR=LL=5 nm (red lines), and LR=LL=7 nm (blue lines) with the same condition in Figs. 2(c) and 2(d). Fx and Fz denote the x and z components of the radiation force, respectively, where Fx has been amplified three times so as to contrast with Fz. From Fig. 2, one can find easily some important properties of resonant radiation force spectra: (i) The resonance peak does not occur at , whereas it has a deviation relative to the energy separation more or less, namely, the resonance peak position occurs at , where is called a radiation shift, which is attributed to the nonlocal optical response in the CQWN system. The magnitude of the radiation shift is dependent on the structure parameters and symmetry of the CQWN. When the CQWN is asymmetric, the radiation shift is not evident, that is, the position of the resonance peak of radiation force occurs at , as shown in Figs. 2(a) and 2(b). When the CQWN is symmetric, however, the radiation shift is evident, and the wider the symmetric CQWN is, the larger the radiation shift will be, as shown in Figs. 2(c) and 2(d). It is because that, compared with asymmetric CQWN, the wave functions of symmetric CQWN would be broadened more due to stronger coupling effects between two wells of symmetric CQWN, which leads the characteristic response scale of electrons to increase. (ii) Fz is evidently larger than Fx in the same polarized state. This is because the quantum-size effect plays an important role in the z direction, which induces the increase of the intensity of electric polarization, due to the charged particles are confined in the z direction of the system.

Fig. 2. Resonant radiation force spectra (Fx dashed lines; Fz solid lines) in the two cases of p- and s-polarization for two different width ratios: LR/LL=2.5/5 (black lines), LR/LL=7.5/5 (red lines) [panels (a) and (b)], and for the three different CQWN widths: LR=LL=3 nm (black lines), LR=LL=5 nm (red lines), and LR=LL=7 nm (blue lines) [panels (c) and (d)]. The mid-barrier width LB=1 nm, and the barrier height V1=V2=2.0 meV, .

In order to contrast the difference of the resonant radiation force under local and nonlocal response, we also plot the resonant radiation force spectra for local and nonlocal cases in p-polarization in Fig. 3. From Fig. 3, we can find that the nonlocality of optical response has an important influence on the resonant radiation force spectra. Obviously, the resonance peak of Fx and Fz has an evident blue-shift under nonlocal response, compared with those under local response with same conditions for p-polarized state. Furthermore, the amplitudes of Fx and Fz under nonlocal response are larger than those under local response, respectively. In addition, for s-polarized state, both Fx and Fz are equal to zero under local response due to the y component of polarization is equal to zero in dipole approximation. However, they are not equal to zero under nonlocal response. The radiation shift of the resonant radiation force spectra due to the nonlocality of optical response is also demonstrated by Eq. (19). From Eq. (19) we can easily see that the resonance structures of the optical force spectra are determined by the roots of det . According to the expressions of the diagonal matrix elements Dxx, Dyy, and Dzz, it is obvious that the resonance peak does not occur at , but has a shift relative to the level interval , which is attributed to the contributions from the real part of the terms , , , respectively.

Fig. 3. Under local (solid lines) and nonlocal (dashed lines) response, the resonant radiation force spectra (Fx and Fz) in p-polarization for the well width: LR=LL=7 nm. The mid-barrier width LB=1 nm, and the barrier height V1=V2=2.0 meV. .

To see how the well-width ratio and the polarized state influence the resonant radiation force spectra, we plot the maximum radiation force (Fx and Fz) at resonance in Fig. 4 and the position of the resonance peak ( ) in Fig. 5 versus the well-width ratio LR/LL in the two cases of p- and s-polarizations, respectively. The black solid, red dashed, and blue dotted lines correspond to the barrier height ratios V1/V2=0.4, 0.7, and 1.0, respectively, with LL=4 nm and V2=2.0 meV. In Figs. 4 and 5, it is evident that both the maximum radiation force and the position of the resonance peak are closely relevant to the structure parameters of the CQWN, and the changing trend of Fx and Fzversus the well width ratio LR/LL is consistent essentially for p-polarization and s-polarization. From Fig. 4, we can see that with the increase of the well width ratio LR/LL, the maximum resonant radiation force decreases at first and then rises; when LR/LL=1, that is, when the CQWN becomes symmetric, it attains a peak value. The CQWN degenerates to a single quantum well when LR/LL=0, and the resonant radiation force is maximal at this point. When LR/ , the resonant radiation force does not increase with further increasing the CQWN width. When the well width ratio value LR/LL is in the approximate range from 1.5 to 1.8, a relative minimal resonant radiation force can be obtained. Furthermore, for a fixed CQWN width ratio, the resonant radiation force becomes smaller and smaller with the increase of the mid-barrier height V1, when the width ratio LR/LL is in the range from 0.5 to 2.0. All the properties manifest that the resonant radiation force can be modulated by designing the structure parameters of CQWN.

Fig. 4. The maximum resonant radiation force (Fx and Fz) versus the well-width ratio LR/LL in the two cases of p- and s-polarizations for V1/V2=0.4 (black solid line), 0.7 (red dashed line), 1.0 (blue dotted line) respectively, with the left well width LL=4 nm, and the outer barrier height V2=2.0 meV, .
Fig. 5. The resonant peak position of optical force spectra, normalized by energy-level interval , versus the well width ratio LR/LL with the V1/V2=0.4 (black solid lines), 0.7 (red dash lines), 1.0 (blue dot lines) in the two cases of p-polarization and s-polarization, respectively, with the left well width LL=4 nm, and the outer barrier height V2=2.0 meV, .

The well width ratio LR/LL is an important factor in influencing the position of resonance peak of the optical force spectra, as shown in Fig. 5. It is worth noting that the position of resonance peak would deviate from the energy interval due to the nonlocal interaction between electrons and laser fields. As the well width ratio LR/LL increases, the deviation of the resonance peak position is more and more evident; when LR/LL=1, the deviation attains to a peak value, and then it decreases gradually. In other words, the radiation shift of optical force at resonance in symmetric CQWN is larger than that in asymmetric CQWN generally. The reason is that, compared with asymmetric CQWN, the wave functions of symmetric CQWN would be broadened more due to stronger coupling effects between two wells of symmetric CQWN, which leads the characteristic response scale of electrons to increase. In addition, the barrier height also has influence on the radiation shift of the position of resonance peak of the optical force spectra in the two cases of p-polarization and s-polarization. It increases with the increase of mid-barrier height V1, and is larger in the p-polarized state.

The change of maximum resonant radiation force (Fx and Fz) with the barrier height ratio V1/V2 is shown in Fig. 6 for three different well width ratios: LR/LL=2.5/5 (black solid line), LR/LL=5/5 (red dash line), LR/LL=7.5/5 (blue dotted line) in the two cases of p-polarization [(a) and (b)] and s-polarization [(c) and (d)]. We can see clearly in Fig. 6 that the maximum resonant radiation forces (Fx and Fz) decrease gradually with the barrier height ratio V1/V2 increasing at first, and then tending to remain stable. It is because that the coupling effect would be weaken with the mid-barrier height increasing. In addition, for a fixed barrier height ratio V1/V2, the change of maximum resonant radiation force (Fx and Fz) with the well width ratio LR/LL is also in agreement with those shown in Fig. 3.

Fig. 6. The change of maximum resonant radiation force (Fx and Fz) versus the barrier height ratio V1/V2 for three different well width ratios: LR/LL=2.5/5, 5/5, 7.5/5 in the two cases of p-polarization [(a) and (b)] and s-polarization [(c) and (d)], respectively, with the left well width LL=5 nm, and the outer barrier height V2=2.0 meV, .

In Fig. 7, we show the maximum resonant force (Fx and Fz) versus the incident angle θ in the two cases of p-polarization [(a) and (b)], and s-polarization [(c) and (d)] for three different well-width ratios, respectively. We note that the maximum resonant radiation forces for the p-polarized state are much larger than those for the s-polarized state, except the case of θ=0. This is because there is a quantum-size effect for the p-polarized state, while it hardly works for the s-polarized state. Furthermore, we can see clearly that the incident angle has an evident influence on Fx and Fz. For p-polarized state, as the incident angle θ increases, it is obvious that both Fx and Fz increase in Figs. 7(a) and 7(b). This is because with θ increasing the z component of the incident field increases, causing an increase of the z component of the polarization in the CQWN, and the quantum-size effect plays an important role in the z direction. When θ =0°, which means that the incident light is perpendicularly incident on the CQWN, and the electronic polarization parallels to the x direction, as a result that Fz at resonance is very small because the quantum-size effect have not any influence in this direction. When θ =90°, that is, the incident light is horizontally incident on the CQWN, the electronic polarization direction is completely parallel to the z direction, which leads to Fx and Fz at resonance to be maximum on account of the quantum-size effects. For the s-polarized state, we can see that decreases and Fx increases with the incident angle θ increasing in Figs. 7(c) and 7(d). The reason for this is that the z component of wave vector will decrease and the x component will increase with the increase of incident angle, respectively. These properties make the CQWN a good candidate for optical manipulation of nano-objects using laser-induced radiation force.

Fig. 7. Maximum resonant radiation force (Fx and Fz) versus the incident angle θ in the two cases of p-polarization [(a) and (b)], and s-polarization [(c) and (d)] for three different well-width ratios, respectively, with the left well width LL=5 nm, and the barrier height V1=V2=2.0 meV, .
4. Conclusion

In this paper, we have theoretically investigated the nonlocal effects on optical force exerted on the AlGaAs/GaAs CQWN from the resonant intersubband transitions for two kinds of polarized states. We have shown that there exists a blue-shift of radiation force spectra at resonance caused by the spatial nonlocality of optical response, and with well width increasing, the radiation shift increases for symmetric CQWN, while it is more evident and larger than those for asymmetric CQWN. Moreover, it is also demonstrated that the resonant radiation force is controlled by the CQWN width ratio, the barrier height ratio, the polarized and incident directions of incident light. The designable CQWN provides a more advantageous way to obtain the suitable optical resonant radiation force. These properties have potential applications in detecting internal quantum properties inside nanostructures, and can open novel and exalting possibilities for the characterization, assembly and optical control of nano-objects.

Reference
[1] Ashkin A 1970 Phys. Rev. Lett. 24 156 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.24.156
[2] Ashkin A 1970 Phys. Rev. Lett. 25 1321 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.25.1321
[3] Mandal S Serey X Erickson D 2010 Nano. Lett. 10 99 https://pubs.acs.org/doi/abs/10.1021/nl9029225
[4] Righini M Ghenuche P Cherukulappurath S Myroshnychenko V García de Abajo F J Quidant R 2009 Nano. Lett. 9 3387 https://pubs.acs.org/doi/abs/10.1021/nl803677x
[5] Nieto-Vesperinas M Sáenz J J Gómez-Medina R Chantada L 2010 Opt. Express 18 11428 https://www.osapublishing.org/oe/abstract.cfm?uri=oe-18-11-11428
[6] Jauffred L Richardson A C Oddershede L B 2008 Nano. Lett. 8 3376 https://pubs.acs.org/doi/abs/10.1021/nl801962f
[7] Chen Y F Serey X Sarkar R Chen P Erickson D 2012 Nano. Lett. 12 1633 https://pubs.acs.org/doi/abs/10.1021/nl204561r
[8] Chaumet P C Nietovesperinas M 2000 Opt. Lett. 25 1065 https://www.osapublishing.org/ol/abstract.cfm?uri=ol-25-15-1065
[9] Maragò O M Jones P H Gucciardi P G Volpe G Ferrari A C 2013 Nat. Nanotechnol. 8 807 https://www.nature.com/articles/nnano.2013.208
[10] Cho S H Lee J Y Han J K Cho B I 2010 Ultrasound. Med. Biol. 36 202 https://www.sciencedirect.com/science/article/pii/S0301562909015786
[11] Li Y Hu Y J 2013 Chin. Phys. B 22 034206 https://cpb.iphy.ac.cn/EN/abstract/abstract52490.shtml
[12] Wang G H Guo Q Wu L J Yang X B 2007 Phys. Rev. B 75 205337 https://journals.aps.org/prb/abstract/10.1103/PhysRevB.75.205337
[13] Wang G H Yan X S Zhang J K 2017 Chin. Phys. B 26 106802 https://cpb.iphy.ac.cn/article/2017/1907/cpb_26_10_106802.html
[14] Wang C Y Wang G H 2014 Chin. Phys. B 23 127103 https://cpb.iphy.ac.cn/EN/abstract/abstract62100.shtml
[15] Liu A Keller O 1995 Phys. Scr. 52 116 http://iopscience.iop.org/article/10.1088/0031-8949/52/1/021/meta
[16] Chang R Leung P T 2006 Phys. Rev. B 73 125438 https://journals.aps.org/prb/abstract/10.1103/PhysRevB.73.125438
[17] Zhang X M Xiao J J Zhang Q 2014 Chin. Phys. B 23 017302 https://cpb.iphy.ac.cn/EN/abstract/abstract57542.shtml
[18] Lu J H Wang G H 2016 Chin. Phys. B 25 117804 https://cpb.iphy.ac.cn/EN/abstract/abstract68531.shtml
[19] Teperik T V Nordlander P Aizpurua J Borisov A G 2013 Phys. Rev. Lett. 110 263901 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.110.263901
[20] Mcmahon J M Gray S K Schatz G C 2009 Phys. Rev. Lett. 103 097403 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.103.097403
[21] Cho K 2003 Optical Response of Nanostructures: Microscopic Nonlocal Theory Berlin Springer-Verlag 1 32
[22] Belov P A Marques R Maslovski S I Nefedov I S Silveirinha M Simovski C R Tretyakov S A 2002 Phys. Rev. 673 181
[23] Mortensen N A Raza S Wubs M Sondergaard T Bozhevolnyi S I 2014 Nat. Commun. 5 3809 https://www.nature.com/articles/ncomms4809
[24] Batke E Heitmann D Tu C W 1986 Phys. Rev. 34 6951
[25] Raza S Bozhevolnyi S I Wubs M Mortensen N A 2015 J. Phys.: Condens. Matter. 27 183204 http://iopscience.iop.org/article/10.1088/0953-8984/27/18/183204/meta
[26] Castaldi G Galdi V Alu A Engheta N 2011 Phys. Rev. Lett. 108 063902 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.108.063902
[27] Wurtz G A Pollard R Hendren W Wiederrecht G P Gosztola D J Podolskiy V A Zayats A V 2011 Nat. Nanotechnol. 6 107 https://www.nature.com/articles/nnano.2010.278
[28] Luo Y Fernande-zdominguez A I Wiener A Maier S A Pendry J B 2013 Phys. Rev. Lett. 111 093901 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.111.093901
[29] Kosionis S G Terzis A F Yannopapas V Paspalakis E 2014 J. Phys. Chem. C 116 23663 https://pubs.acs.org/doi/abs/10.1021/jp3090183
[30] Jackson J D 1999 Classical Electrodynamics New York Willey 3